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PREFACE 


This  research  report  was  prepared  by  Frank  C.  Liu,  Professor  of 
Mechanical  Engineering,  The  University  of  Alabama  In  Huntsville  during  his 
leave  at  the  Air  Force  Institute  of  Technology.  The  objective  of  this 
research  is  to  Investigate  dynamic  response  of  a  tether  connected  satellite 
system  due  to  mass  transport  along  the  tether.  This  problem  may  have  use¬ 
ful  application  to  NASA  "Skyhook"  program.  The  author  is  regretful  that 
due  to  limited  time  and  computer  programming  skills,  and  the  large  number 
of  computer  runs  required,  that  numerical  results  are  not  available  at  the 
present  time.  They  will  be  presented  in  a  later  report. 
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ABSTRACT 


^  Two  satellites  connected  by  a  long  flexible  tether  along  the  earth 
radial  direction  comprise  a  stable  equilibrium  state.  This  research 
investigates  the  manner  in  which  a  third  mass  transporting  from  one 
satellite  to  the  other  disturbs  the  equilibrium  state.  A  system  of  four 
equations  of  in-plane  motion  has  been  derived  based  on  the  assumptions  that 
the  tether  remains  straight  between  the  masses  and  of  constant  length.  A 
combination  of  computer  subroutines  DGEAR  and  ZANAYT  is  suggested  for  the 
approximate  solution  of  the  system  of  nonlinear  differential  equations  with 
one  constraint  condition  on  the  variables.  Alternatively,  a  system  of 
four  independent  differential  equations  are  derived  by  eliminating  the 
Lagrange  multiplier,  n 


PREFACE 


ACKNOWLEDGEMENTS 

ABSTRACT 

TABLE  OF  CONTENTS 
NOMENCLATURE 

I.  INTRODUCTION 

II.  ANALYSIS 

1.  Formulation  of  Kinetic  and  Potential  Energies  of  the 
Satellite  System 

1.1  Coordinate  and  Definition  of  Variables 

1.2  Formulation  of  the  Kinetic  Energy  of  the  Masses 

1.3  Formulation  of  the  Potential  Energy  of  the  Masses 

1.4  Formulation  of  the  Kinetic  Energy  of  Tether 

1.5  Formulation  of  the  Potential  Energy  of  Tether 

2.  Lagrange's  Equations  of  the  System 

2.1  Equation  of  Constraint 

2.2  Lagrange's  Equation  of  Motion 

2.3  Equations  of  Motion 

3.  Generalized  Applied  Forces  and  Initial  Conditions 

3.1  Free  Motion  of 

3.2  Forced  Motion  and  Generalized  Forces 

3.3  Reversed  Position  of  the  Satellite  System 

3.4  Initial  Value  of  x* 

4.  Numerical  Methods  of  Solution 

4.1  Integration  by  Combined  Subroutines 

4.2  Elimination  of  the  Lagrange  Multiplier 

4.3  An  Example  of  Illustration 


4.4  Numerical  Verification 


32 


4.5  Parametric  Study  34 

4.6  Computer  Print-out  and  Computer  Time  Control  34 

III.  RECOMMENDATIONS  AND  CONCLUSION  35 

IV.  APPENDICES  37 

A.  Definition  of  Notation  Used  in  Tables  1  and  2  37 

B.  Definition  of  Coefficients  of  Equations  of  Motion  38 

C.  Computer  Subroutines  41 

V.  REFERENCES  42 


v 


NOMENCLATURE 
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in  Apendix  B 
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defined  by  Equation  (4.3) 

defined  by  Equation  (4.5) 

constraint  equation.  Equation  (4.2e) 

defined  in  Appendix  A 

unit  vector  along  orbital  radius 

unit  vector  along  orbital  velocity 

length  of  tether 

=  m^  +  +  m^ 

=  m/m^ 

defined  by  Equation  (1.5) 
mass  of  outer  satellite 
mass  of  inner  satellite 
mass  of  transport  mass 
mass  of  tether 
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NOMENCLATURE 


defined  by  Equation  (1.3) 
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2 

=  p(t)/m3tu»0 

defined  by  Equation  (3.13) 

min.  magnitude  of  driving  force  to  start  transfer  motion 

generalized  coordinate 

generalized  force  corresponding  to  q. 

defined  by  Equations  (2.8),  (3.7)  and  (3.13) 

subscript  designates  the  corrrespondi ng  equation, 

defined  by  (4.3) 

defined  by  Equation  (4.4) 
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position  vector  of  from  center  of  the  Earth, 
i  =  1,  2,  3 

position  vector  of  center  of  mass  of  satellite  system 
*  r0/t 

integration  variable  along  tether 
time  variable 

kinetic  energy  terms  defined  in  Table  1. 
kinetic  energy  of  masses 
kinetic  energy  of  tether 


Va  t0  Vi 


potential  energy  terms  defined  in  Table  2. 

=  -  u/r.j 

velocity  vector  of  m^ ,  i  =  1,  2,  3 
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potential  energy  of  masses 
orbital  velocity  of  center  of  mass 

=  -  u/rQ 

potential  energy  of  tether 

coordinates  defined  in  section  4.4 

variable  of  equations  of  motion  defined  i 
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defined  by  Equation  (1.3) 

=  a  -  y,  angle  between  p  and  p« 
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=p  2n 

coordinates  defined  in  section  4.4 

angle  between  p ^  and  rQ 
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2 

=  x/nuui  ;  also  used  in  Section  4.3 
o 

gravity  constant  of  Earth 

position  of  vector  of  m. ,  i  a  1,  2,  3 


=  P/t 


=  uiQt,  non-dimensional  independent  time  variable 
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SYMBOLS: 


a  letter  underlined  denotes  vector 

dot  between  vectors  denotes  dot  product 

dot  on  top  of  a  letter  denotes  time  derivative 

double  dots  on  top  of  a  letter  denotes  second  time 

derivative 

prime  denotes  derivative  with  respect  to  x 

double  prime  denotes  second  derivative  with  respect  to 

normalized  quantity 

summation  of  all  variables  in  the  equation 


I.  INTRODUCTION 


Since  Colombo  ^  developed  the  concept  of  connecting  a  heavy  mass  to 
a  satellite  by  a  very  long  flexible  tether  in  1974,  the  subject  has  stimu¬ 
lated  the  science  and  engineering  community.  A  recent  article  by  Ivan 
Bekey  of  the  NASA  Office  of  Space  Flight  gives  detailed  descriptions  of 
many  scientific  applications  of  this  idea.  Two  contractors,  Martin 
Marietta  (Denver),  Aerospace  and  Ball  Aerospace,  have  been  given  respon¬ 
sibility  for  the  design  of  the  "skyhook"  project  which  is  scheduled  for 
first  flight  by  NASA  in  1987. 

Preliminary  analyses,  feasibility  studies,  and  design  of  a  Tethered 
Satellite  System  (TSS)  were  conducted  at  Marshall  Space  Flight  Center  by 
Rupp  and  Lane  and  Baker,  et  al  Numerous  papers  have  been 

published  in  the  last  decade  dealing  with  the  dynamics  of  deployment  and 
retrieving  the  mass  from  a  space  shuttle.  Various  dynamical  models  have 
been  developed  by  many  investigators.  These  models  may  include  one  or  more 
of  the  following  actions: 

(a)  tether  mass, 

(b)  three-dimensional  motion, 

(c)  longitudinal  vibration  of  the  tether, 

(d)  transverse  vibration  of  the  tether, 

(e)  rotational  motion  of  masses, 

(f)  offset  distance  of  attachment  to  C.M.  of  masses,  and 

(g)  eccentricity  of  TSS  orbit. 

r  5  n 

Comparisons  of  the  models  can  be  found  in  the  paper  by  Misra  and  Modi  L 
A  formulation  of  a  general  dynamical  model  for  TSS  by  the  same  authors  is 
given  in  Reference  6. 


Another  aspect  of  the  problem  Is  optimal  control  of  the  tension  in  the 
tether  for  dynamic  stability  daring  deployment  and  retrieving  of  the  mass. 
See,  for  example,  a  paper  by  Bainum  and  Kumar 

The  objective  of  this  research  concerns  with  a  problem  which  is  quite 
different  from  that  described  above.  The  TSS  aligned  along  earth  radial 
direction  is  a  stable  equilibrium  state.  Consider  the  requirement  that  a 
third  mass  must  be  transported  from  one  sub-satellite  to  the  other  along 
the  tether.  The  mass  transfer  operation  can  be  accomplished  by  free  motion 
with  sufficient  initial  velocity,  or  by  applying  a  small  thrust.  It  is 
desired  to  find  the  dynamic  response  of  the  TSS  due  to  the  motion  of  the 
transport  mass. 

The  dynamical  model  to  be  treated  here  will  include  only  one  of  the 
factors,  (a),  mentioned  above.  Hence,  the  TSS  has  three  degrees-of-freedom 
There  is  no  convenient  method  for  reducing  the  dynamical  model  to  three 
Independent  variables.  Thus,  one  constraint  condition  must  be  induced  bet¬ 
ween  the  four  variables  used  for  the  formulation  of  equations  of  motion. 
This  constraint  condition  is  that  the  tether  remains  straight  between 
masses  and  of  constant  length.  Due  to  inclusion  of  the  mass  of  the  tether, 
the  system  of  four  second  order  differential  equations  will  be  very 
lengthy. 

To  the  author's  knowledge,  there  is  no  existing  computer  subroutine 
for  the  approximate  solution  of  a  system  of  differential  equations  with  a 
constraint.  A  combined  computer  subroutine  for  the  solution  of  such 
problems  Is  suggested,  and  will  be  tested  on  an  example  with  known  solu¬ 
tion.  Further,  a  new  system  of  four  independent  second  order  differential 
equations  has  been  derived  from  the  original  system  by  eliminating  Lagrange 


II.  ANALYSIS 


I.  FORMULATION  OF  KINETIC  AND  POTENTIAL  ENERGIES  OF  THE  SATELLITE  SYSTEM 
The  following  assumptions  have  been  made  for  the  formulation  of  the 
kinetic  and  potential  energies  of  the  three  masses  and  of  the  tether  con¬ 
necting  them. 

(a)  The  motion  of  the  center  of  mass  of  the  system  is  undisturbed 
by  the  relative  motions  of  the  masses  and  the  tether,  i.e.,  the  center  of 
mass  maintains  uniform  motion  in  a  circular  orbit. 

(b)  The  tether  is  flexible  and  inextensible.  It  remains  straight 
at  all  times  between  the  transport  mass  and  each  of  the  others. 

(c)  The  transport  mass,  nij,  is  free  to  move  along  the  tether, 
i.e.,  friction  is  neglected. 

(d)  There  is  no  motion  normal  to  the  orbital  plane. 

(e)  The  main  satellite  and  the  sub-satellite  m^  are  treated  as 
point  masses,  I.e.,  the  off-set  distances  from  the  attachment  point  to 
center  of  masses  are  neglected  in  the  development  of  the  equations  of 
motion.  However,  satellite  dimension  will  be  partially  accounted  for 
through  the  distance  a  (See  Fig.  1)  in  the  initial  and  final  radii. 

(f)  For  the  same  reason,  the  rotational  inertias  of  all  masses 
are  neglected. 

(h)  Mass  1*1  is  much  greater  than  m^  and  both  mass  m3  and  the  mass 
of  the  tether  m^  are  much  smaller  than  m^. 

(i)  The  driving  force  applied  on  m^  had  negligible  effect  on  the 
orbital  motion  of  the  system. 

1.1  Coordinates  and  Definition  of  Variables 

Figure  1  illustrates  a  tether  connected  satellite  system.  The 
left  side  figure  shows  a  disturbance  from  the  original  undisturbed  con¬ 
figuration  (right).  A  rotating  coordinate  system  is  chosen,  with  origin  at 


the  center  of  mass  of  the  system,  S-axis  along  the  radial  direction,  and 
n-axis  in  the  direction  of  the  orbital  velocity.  Let  be  the  position 
vectors  of  m^  (i  =  1,  2,  3)  relative  to  the  rotating  coordinate  system. 
Denote  by  p  the  position  vector  of  m^  relative  to  m^  and  by  8  and  y  the 
angles  of  p2  and  p  with  respect  to  5 -axis. 

The  tethered  satellite  has  three  degrees-of -freedom.  Two  sets  of 
variables,  (  p2>  9)  and  (  p,y)  are  chosen  for  formulating  equations  of 
motion  of  the  system.  Hence,  the  four  variables  must  satisfy  the  assumed 
constraint  condition  that  the  length  of  tether  is  constant.  Now,  set  the 
vectors 

-2  =  p2  (-cos9i  +  sin0D  (1.1a) 

g  =  p  (cosyJ_-  sinYj_)  =  03  -  P2  (1.1b) 

and  express  p^  and  g3  in  terms  of  g2  and  p. 

It  follows  from  assumption  (a)  that 

mlpl  +  ^  02  +  m3-3  +  mtp  ^-2  *  03^/24  +  m^U  -  p ) ( p ^  +  g3)/24  =  0 
This  equation  gives 

-  Pi  =  i[M23  +  ^Mt{1  +  a)  02  +  (M3  +  ^Mt)pl  [1  +  1/2Mt(1  ’  a)]"1  (1- 

where 

m  =  m^  +  m^  +  m^,  =  nu /m.^  (i  =  1 ,2 ,  t ) 

M  =  m/m1 ,  R23  =  M2  +  M3  =  M  -  1  a  =  p/4 

After  eliminating  second  and  higher  order  terms  of  Mt,  one  may  write 

P1  =  "  ^23-2  +  ^3  (1.3) 

where 

fi23  ’"23(1  ♦»23>  *3  ‘  M3(I  *  “31 

®23  *  Mt(Mo  *  1  -  M23,/2H23  a3  *  V*f  *  1  '  V/2M3 

Note  that  fl23  and  fl3  are  variables  through  inclusion  of  p/4  =  o  in  the 

expressions  for  «23  and  a3.  Due  to  assumption  (h),  the  quantities  which 
involve  Mt  as  much  smaller  than  unity. 


1.2  Formulation  of  the  Kinetic  Energy  of  the  Masses 
The  velocity  vector  of  mass  m^  is 

vt  =  vo  +  £j  +  ^0  x  £,  1  3  1.  2,  3  (1.4 

Substituting  Equation  (1.1)  and  (1.2),  results  in 

Vj  =  ^3  £p2cose  -  p2(9  -  u>0 ) si ne ]  -  [M3  pcosy  -  p ( y  -  u>Q ) si rry ] 

+  {rQu)o  -  fi23[p2sine  +  p2  (0  -  a>o)cos0]  +  M-^psirry  +  p(y  -  u)o)cosY]}j_ 

V2  =  C p2( 9  -  u>o)sin0  -  p2cos0]J_  +  [roa>o  +  p2sin0  +  p2(0  -  ujo)cos0]j_ 

v^  =  Vg  +  Cpcosy  -  p(y  -  (*Jo)sinY]_i_  -  [psirry  +  p(y  -  u>Q)cosY]j_ 

It  is  helpful  to  present  in  tabulated  form  as  in  Table  1  the  individual 
terms  which  constitute  the  kinetic  energy. 

More  mass  parameters  are  defined  in  the  following  for  the  for¬ 
mulation  of  the  kinetic  energy  terms  in  Table  1. 

for  Ta:  Mfl  *  +  m2  +  m3)/m  *  M23(l  +  2a23M23/M)  (1.5a] 


for  T.  :  M.  =  (m.R.,  +  m,)/m,  =  1  +  M,  +  2a0M, 
00  lood  3  33 


(1.5b] 


for  Tc:  Mc  =  (m1R23M3+m3)/m  =  M3  [1  +  (a3  +a23)M2;J/M]  (1.5c] 

It  is  helpful  to  list  derivatives  of  the  above  mass  parameters  which  will 
be  used  in  formulation  of  Lagrange's  equations. 
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*  Mt(l  -  l/2M)M3/t 

Mc  = 

Mt(l  -  1/2M)M35 

(1.5f) 

From  Table  1  the  various  velocity  terms  can  be  formed  by  summing 
up  the  products  of  the  elements  in  the  corresponding  column  and  that  in  the 
first  column.  For  example. 
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The  total  kinetic  energy  of  the  three  masses  is 


2T  =  ntf  T  +  m,MKT.  -  nM  (T  ♦  T.) 
m  aa  3bb  cc  a 


mooao-iT  +  m-a^T,  +  mT  (1.6) 
23  23  e  3  3  f  0 


In  Table  1  the  second  column  from  the  right  is  2T  and  the  last  column 

m 


gives  all  terms  in  the  kinetic  energy  of  the  tether  which  will  be  for¬ 
mulated  later. 

1.3  Formulation  of  the  Potential  Energy  of  the  Masses 

Denote  V.  *  -  u/r^ ,  where  r^  is  the  distance  from  center  of  the 
earth  to  m^  and  u  is  the  gravitational  constant.  The  potential  energy  of 
masses  is  then 

Vm  =  ffllVl  +  Vz  +  Vs 

Let 


(1.7) 


l  =  [( r  +  p,).(r  ♦  p,.)]-1^ 


i 


-o  -1 


-o  -1 


(1.8) 


Making  use  of  Equations  (1.1),  (1.2)  and  the  relationships, 
Co-62*  ‘  ro  p2cos9  * 


r  .0=  rQ  ccosy 


r1  =  ro  Cl  *  (2Co-2i  *  °?,/ro1 


and  taking  the  binomial  expansion. 


i  - 1  Cl 

r1  ro 


(2r  .p,  *p2)/2r2  *  3<r „.P,/r?)/2] 


'-oJ£)l  '  Hi "  '  o  ,-o’=r  '  o' 


one  obtains  the  following  expressions 


j-  =  £  {1  +(^3pcosy  -  fl23  p2cos0)/rQ  +  (3cosZ0  -  l)/2r 


1 


2 

0 


2  2  2  2  2 
+  fi3p  (3cos  y  -  l)/2rQ  -  fl3fl23P2p[  cose  cosy  +  cos(e+  Y)]/rQ}  (1.8a) 


\  ■  p  {  1  +  ■;  P2cose  +  P2(3cos29  -  l)/2r2} 


(1.8b) 


-  s  r  i1  +  r  (p2cos9  '  pC0SY'  +  p2(3cos  9  “  l’/2r0 
3  o  o 

2  2  2  2 
+  p  (3C0S  Y  -l)/2r  -  P2pEc0S9C0SY  +  C0S(9  +  Y )]/rQ} 


(1.8c) 


Individual  terms  constituting  potential  energy  are  tabulated  in 
Table  2.  Thus,  the  potential  energy  of  the  three  masses  can  be  written  in 


the  form 


v™  =  ”Vo  -  m23“23Va  *  "3“3Vb  *  ""a,3Vc  ‘  V 

*  "'3V3ve  ■  V  -  mVvg  *  V 


(1.10) 


where  the  V's  in  Equation  (1.10)  are  defined  in  Table  2.  The  last  column 
of  Table  2  gives  terms  in  the  potential  energy  of  the  tether  which  will  be 
formulated  later. 

1.4  Formulation  of  Kinetic  Energy  of  the  Tether 

It  is  assumed  that  the  tether  remains  straight  between  the  masses 
as  shown  in  Figure  2.  Thus,  velocity  of  points  Pj  and  can  be  written  in 


the  forms. 


vp  =*  [s^  +  U  -  p-  s^Vjl/U  -p)  O^s^t-p  (1.11a) 


Vp2  =  [s2v3  +  (p  -  s2)v2]/p 


0  s  s2  <  p 


(1.11b) 


The  kinetic  energy  of  the  tether  can  be  obtained  by  integration 


2T  =  (mf/t)  /*"p(v  .v  )ds,  +  /p(v  .v  )ds? 

1  t  0  "P1  “P1  1  0  “p2  ~p2  2 

2  2  2 

=  (mt/3i)[(t  -  p)v1  +  pv2  +  iv3  +  (t-pl^.Vj  +  py2.y3]  (1*12) 

Making  use  of  Table  1,  Equation  (1.12)  yields 

2Tt  =  (mt/6)  [E1Ta  ♦  +  G3(Tc  +  T,)  ♦  3(1  -  M23  ♦  Ma)Te 


-  3(1  -  M,  +  M,o)Tr]  +  m.T 


(1.13) 


(  -m  6r  )p?Pcos (e -y ) 


where  the  coefficients  are  defined  in  Appendix  A.  Note  that  the  bar  on  top 
of  M,  M^,  and  in  Equation  (1.13)  and  in  the  expressions  for  its  coef¬ 
ficients  have  been  removed  to  retain  only  first  order  terms  of  This  is 
justified  through  assumption  (h). 

1.5  Formulation  of  the  Potential  Energy  of  the  Tether 

The  potential  energy  of  the  tether  may  be  written  in  the  form 


Vt  =  (um./A)  (  /*~p(l/r  Ids.  +  /p  ( 1/r  )ds?| 

x  z  0  P1  1  0  Pz 


(1.14) 


where  r  and  r  are  respectively  distances  of  P,  and  P„  from  the  center 
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of  the  earth,  as  shown  in  Figure  2.  The  position  vector  of  points  P^  and 
?2  on  the  tether  may  be  expressed  in  the  forms 

rp  =  [s1r1  +  U  -  ps^r-jj/U-p) 


-p2  =  Cs2-3  +  (p‘  s2J-2^/p 
l/r_  =  ( r_  .r_  )’^2  i  =  1,  2 


Pi  -Pi  -P1 

Using  the  following  relationships. 


(1.14a) 

(1.14b) 

(1.14c) 


2  2  2 
r-'  =  ro  +  2-o*^i  +  Pi 


1 


one  obtains 


=  ro +  c0-fii  +  Vfij  +  fii-fij 
?  =  J  l1  -  (r0*ei)[si  +  si  u  •  p_  si)]  '  (r0-e3)  x 

pi  0 

2  2  2  2 
C(t  -  p-  Sj  )  +  s1  ( 2  -  p-  sx )]  +  3(rQ.p^)  [sj  +  s^t  -  p-  Sy  )/2rQ 

2  2  2 
+  3(rQ.pj)  [(2  -p  -  Sl)  +  s1  (2  -  p-  S|_ ) ]/2rQ 

2  2  2  2  2  2 
-  Pj  +  (A  -  p-  Sj  )  P3  +  Sj  (t  -P-  Sj  )g1  .P3]}/r0  (2  -p) 
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^  ^  i  1  "  *  s 2 ^  P—  S£ ) ]  "  ( f 0 •  )  C(  P-  ^2  )  ■*■  $2(P-  ^9^ 

p2  0  "  "  '■'  ^  c 


2  2 


+  3( rQ. p3)  [s2  +  s2(p  -  s2)]/2rQ  +  3(r0.o2)  [(P  -  s2)  +  s?(P  -  s2)]/2r‘ 


22  22  22 

-  1/2  [Sg  p3  +  (p  -  s2  )  P2  +  s2  (p  -  ^2  )p2  .£3  ]  It  p  } 


The  integrals  lead  to 


l‘-j  1S1  >rl  ds2  -  1  11  -  1  [*Pj  ♦  PP2  P  P(P2.P3)  *  U  -P)(P[.P; 
0  Pi  0  P2  0  6ro 


• 2  L  -0  -3 


-o-l 


1  r  2  2  2 

+  Wr0*V  +  p(r0-p2)  +  ^  •  p)(r0.pi  )  ] 

2ro 

+  _2  ^-0*23^  ^-0*22^  +  ^  -p'(  W1} 

ro 

Substitution  of  the  following  into  the  above  equation 
2  22  22 

P1  =  ^23  p2  +  M3  p  '  2M3M23P2Pcos  (0  "  Y^ 

222 

P3  =  P2  +  p  -  2p2pcos(9  -  y) 

-o*2l  =  VM23p2Cos9  "  M3  PC0Sy^ 

r  . p0  =  -  r  p0cos9 
-0  -Z  02 

rn. 0,  -  r  (pcosy  -  p„cos0 ) 

-0  -J  0  c 

results  in  the  terms  given  in  the  last  column  of  Table  2.  Thus,  the  poten 
tial  energy  of  the  tether  is 

Vt  =  5mt  (3(1  “  M23  +  ‘b)Va  -  '  M3  +  M3°)vb  +  3Eivc  *  2(1  +  2o)Vd 


+  3GlVe  -  2Vf  +  E2Vi>  +  Vo 


(1.15 


where  the  Vs  are  defined  in  Table  2. 


2.  LAGRANGE'S  EQUATIONS  OF  THE  SYSTEM 


It  is  straightforward  to  apply  Lagrange's  method  to  obtain  equations 
of  motion  in  terms  of  the  chosen  variables  (p2>  9)  and  (p,y).  The 
constraint  that  the  length  of  the  tether  must  be  a  constant  may  be 
satisfied  by  introducing  a  Lagrange  multiplier  into  the  formulation. 

2.1  Equation  of  Constraint 

From  Figure  1,  it  can  be  seen  that 

I  Pj-  P3  I  +  I  P  I  =  i  (2.1) 

Making  use  of  Equation  (1.3),  Equation  (2.1)  becomes 

[(1  +  f^3)  p2+  (1  +  M3 ) p] . [ ( 1  +  ^23)P2  +  (1  +  M3)P]  =  U  -  P)2  (2.2) 

The  above  equation  yields 

fl%2  +  (1  +  M3)2p2  -  2R(1  +  R3)p2Pcos(0  -  y)  -  U  -  P)2  =0  (2.3) 

Direct  differentiation  of  Equation  (2.3)  leads  to 

a  p,  +  ano  +  a.H  a  y  =  0  (2.4) 

P2  2  p  9  y 

where,  after  second  and  higher  order  terms  of  Mt  have  been  eliminated, 


a 

P2 


[Mp2  -  (1  +  M3)  pcos (9  -  y)  +  Mt^aap2  “  0cp)]/4 


ap  =  +  ~  Up/M  +  i/M  -  (1  +  M3)p2cos(9  -  y) 

+  (Mt/2t)[Mp2  +  (1  +  M3)(M3/M)p2  -  (1  +  2M3)PP2 
+  2*{abP  -  acP2  )]}  /ft 
aQ  =  (1  +  M3  +  Mtac)P2osin(0  -  y) 


(2.5a) 


(2.5b) 

(2.5c) 


a..  =  -  (1  +  M-,  +  M,.a  )P„asin(9  -  y) 


(2.5d) 


In  the  above  expressions  the  following  notations  have  been  used 


aa  =  Mo  +  1  - 

ab  =  (M3a  +  1  ‘  M3)(1  +  M3)/M 
ac  =  (Vfc  +  Hj)a  +  (1  -  HjM23  )/M 

Note  that  Equation  (2.4)  has  been  divided  by  a  factor  2Mt  to  obtain 
Equation  (2.5)  which  will  give  the  Lagrange  multiplier  used  in  the 
equations  of  motion  a  dimension  of  force. 

2.2  Lagrange's  Equations  of  Motion 

Let  the  kinetic  energy  and  potential  energy  of  the  satellite 
system  respectively  be 


T  =  Tm  +  %  , 
m  t 


V  =  Vm  +  V* 
m  t 


where  the  terms  on  the  right  hand  side  are  given  by  Equations  (1.6),  (1.7), 
(1.12)  and  (1.15).  Four  equations  of  motion  of  the  satellite  system  are 
obtained  from  Lagrange  equation, 


J  iL 

dt  3^ 


3T  t  3 Y 
3q.  dq. 


(2.6) 


q.j  -  P2 »  P »  0  »  and  Y 

where  a^  are  given  by  Equation  (2.5)  and  are  generalized  forces  due  to 
force  applied  on  m^.  Lagrange  multiplier  X  is  introduced  as  a  result  of 
the  constraint  condition  given  by  Equation  (2.4). 

2.2.1  Elimination  of  Nonlinear  Terms 

Note  that  the  variable  p  which  represents  the  distance  of 
m^  from  mg  varies  from  i  to  zero  during  the  transfer  and  is  the  only 
variable  which  cannot  be  treated  as  a  first  order  small  quantity.  Usually 
the  variable  pg  is  replaced  by  its  initial  equilibrium  length  plus  a 
displacement  variable;  but  since  the  system  equations  of  motion  cannot  be 
linearized  due  to  p,  will  also  be  kept  as  a  finite  variable.  The  deri- 


sidered  as  first  order  small  quantities.  All  the  second  and  higher  order 
quantities  will  be  neglected  from  the  equations  of  motion. 

2.2.2  Normalization  of  Variables 

All  the  mass  quantities,  m,  n^,  m^,  and  m^  are  normalized  in 
terms  of  the  outer  satellite  m^  and  are  denoted  by  their  capital  letters. 
All  the  length  variables  are  normalized  in  terms  of  the  length  of  the 
tether,  z,  as  follows: 

P2/i  =  <S  ,  p/z  =  o,  and  rji  =  Rq  (2.6) 

The  time  derivatives  of  variables  are  normalized  by  the  orbital  period. 
Setting 

3  1 

(p/rQ)  ,c  =  a)  and  t  =  u)Qt 

one  may  change  the  time  derivatives 

-  2 

$i  =  ^oqJ  q1  *  wQq?  (q.  =  P2,  P,  9.  and  y)  (2.7) 

where  the  prime  denotes  differentiation  with  respect  to  t. 

To  perform  the  normalization  process,  the  four  equations 

obtained  by  the  Lagrange's  method  for  the  variables  P2,  P,  9,  and  y  are 

2  2  2  2 

divided  respectively  by  m^JUQ,  m^Zu^,  mitu,0P2*  and  ntjZu)op.  Thus,  all  the 
variables  and  their  coefficients  in  the  equations  of  motion  are  dimen¬ 
sionless. 

2.2.3  Change  Variable 

It  is  more  meaningful  to  use  6  =  9  -  y  to  replace  y  as  a 
dependent  variable,  while  the  angle  y  is  used  only  for  the  convenience  of 
forming  the  energies.  As  shown  in  Figure  1,  8  represents  the  angle  between 
the  vectorsp2  and  p.  The  initial  and  final  values  of  8  are  always  zero. 


2.2.4  Remark  on  the  Differentiation  of  Enery  Terms 

It  is  important  to  note  that  all  the  coefficients  of  the 
energy  terms  given  In  Tables  1  and  2  vary  with  a.  Therefore,  in  for¬ 
mulating  Lagrange  equations  differentiations  must  be  carried  out  on  energy 
terms  as  well  as  their  coefficients.  However,  in  the  coefficient  of  the 
energy  terms  for  the  tether  all  variable  mass  parameters  ft,  f^,  and  ft^  are 
replaced  by  their  constant  counterpart  M,  M23,  and  M3  respectively,  since  all 
tether  energies  already  have  a  first  order  coefficient  M^.  See  assumption 
(b). 

2.3  Equations  of  Motion 

Four  equations  of  motion  obtained  from  Lagrange's  equations  are 
presented  in  the  forms 


Cu6>  .  Cl2o" 

-  [A^5  + 

A13°  *  V  + 

‘IS*'1  =  V*  * 

(2.8a) 

C21S"  *  C22°" 

- c  V 

+  a23.  .  A26e- 

+  A28s']  =  A29A*  ♦  QJ 

(2.8b) 

°119"  *  °12  6* 

-  [8126' 

*  ei4«'  ♦  B158 

*  B178]  ■  Bj9A*  ♦  Q*6 

(2.8c) 

V  *  °226‘ 

-  [8225' 

*  824<j'  *  B258 

+  B27al  ’  B29A*  *  °*3 

(2.8d) 

where 

x*  =  x/m  Q*2  =  QpjM<V  Q*  =  Q9M^p2, 

%  =  Wv  and  QV  VVv 

Next,  rewriting  the  equation  of  constraint  given  by  Equations  (2.3)  and 
(2.4)  In  terms  of  the  non-dimensional  variables  gives 

ftV+  (I  +  ft3)V  -  2ft  (I  +  M3)  o6cos3  -  (1  -  o)2  =  0  (2.8e) 

•  2 

[M«  -  (1  +  M,)oC0S3  +  M*  (a  6  -  a  a)]  5  +  I  [ ( 1  +  M- )  -  l]o  +  1 


and 


-  M(1  ♦  M3)6COS0  +  l/2Mt CM  5  +  (1  +  MjJMjCf  -  M(  1  +  2Kj  )<j<5 


+  2M(a.  a  -  a  6)]  a/M  +  3  M(  1  +  M-  +  Mfca  )o5sin6  =  0  (2.8f) 

DC  o  t  C 

The  solution  of  the  four  equations  of  motion  must  satisfy  the  constraint 

equation  given  by  either  (2.8e)  or  (2.8f). 

Note  that  at  the  initial  equilibrium  state,  6  *  8"  *  S'  =  o'  =0,  and 

Equations  (2.8c)  and  (2.8d)  both  reduce  to  a  single  equation. 

9"  +  3a  =  0  (2.9) 

[81 

This  equation  represents  oscillations  of  a  dumbbell  in  a  circular  orbit1  . 
2.3.1  Equations  of  Motion  With  =  0 

The  system  of  equations  of  motion  can  be  simplifed  con¬ 
siderably  if  Mt  is  set  to  zero.  Equations  (2.8)  reduces  to 

M236"  "  M3°"  +  2(M236  ■  M3o)tt'  +  2M3a3'  -  3M235  +  3M3a 


=  X*[M<5  -  (1  +  M-)oCOSB]  +  Q* 

J  Pr 


(2.10a) 


(1  +  M3)o"  -  M6"  +  2  [(1  +  M3)o  -  M6]0'  -  2(1  +  M3)aB '  +  3(1  +  M3)a  +  3M6 

=  X*  {[(1  +  M3)2  -  l]a  +  1  -  M( 1  +  M3)6cosb}  M/M3  +0*  (2.10b) 

(M236  -  M3o)(6-  +  39)  +  M3oB"  -  2M236‘  +  2M3a‘ 

=  X*( 1  +  M3)asinB  +  Q*q  (2.10c) 

[(1  +  M3)a  -  M6 ] ( 9 *  +  39  -  36  )  -  (1  +  M3)o6"  +  2M6 '  -  2(1  +  M^o' 

=  -  X*  [(1  +  M3)M/M3]  5 sine  +  Q*&  (2 . lOd) 


It  is  important  to  find  out  numerically  whether  neglecting  the  contribution 
of  the  tether  mass  leads  to  any  significant  differences. 


3.  GENERALIZED  APPLIED  FORCES  AND  INITIAL  CONDITIONS 


Two  classes  of  mass  transfers  are  considered  in  this  paper.  One  class 
of  transfer  is  that  m^  is  given  an  initial  velocity  with  sufficient  magni¬ 
tude  to  cross  ' ne  orbit.  In  this  case  all  the  generalized  forces  in 
Equation  (2.8)  are  zero.  The  other  class  of  mass  transfer  is  that  m^  is 
driven  by  a  continuous  thrust  with  sufficient  initial  magnitude  to  get 
motion  of  m^  started.  The  free  motion  of  the  transfer  mass  is  treated 
first. 


3.1  Free  Motion  of 

Motion  in  this  type  of  transfer  depends  entirely  on  initial  con¬ 
ditions.  Two  cases  are  considered,  one  is  that  m^  departs  from  the  outer 
satellite,  m^,  and  the  other  is  that  m^  starts  from  the  inner  satellite,  m2 
3.1.1  Departure  From  m^ 

To  derive  the  necessarily  initial  value  of  6q,  disregard  the 
off-set  distance  from  center  of  m^  at  starting  point,  From  definition 
of  center  of  mass  of  the  system, 

1/2  m^z  +  (m^  +  iUjH  =  (m  +  mt)  p 2 ( 0 ) 


It  gives 


5o  ■  11  ♦  "3  *  ^  V'V  Mmt  ’  "  +  Mt 

Now,  let 


(3.1a 


o0  «  1  -  Uj/t)  (3.1b) 
The  minimum  Initial  velocity  of  m^  to  cross  the  orbit  is  determined  by  work 
and  energy  principle.  It  follows 


,  2  ,rl  2  2  ,2  22 

V2  =  /  (ru»0  -  u/r  )dr  =  l/?^  -  rQ  hQ  +  y/^  -  \x/rQ 

ro 


3  ^2  rn  C(  1  ■  Pon  ! rn  ) 


♦ 


[d  + 


P?n 


_i 


-  Hu/r„ 


Substituting  Equation  (3.1a)  and  changing  to  dimensionless  form,  it  gives 

o;  <  -  V?  [1  -  (1  ♦  iy  \  Mt)/Mmt]  =  -  ,’3  (M2  .  \  «t)/«mt  (3.3) 

The  negative  sign  shows  that  the  velocity  is  in  the  direction  of  decreasing 
a  as  m^  moves  toward  i^. 

3.1.2  Departure  From  m^ 

The  initial  position  of  the  three  masses  gives 
1/2  mt*  +  =  (m  +  mt)P2(0) 


Hence, 


so  -  11  ♦  i  V/Mmt  «o-V‘ 


(3.4) 


Similarly,  applying  the  work  and  energy  principle, 

1  .2  i  2  2  2 

2  po*  2  r2”  'o^o  +  p/r2'w/ro 

=  1/2  rQJ0Ul  -  P20/ro)2  -  1]  ♦  [(1  -p2Q  /r0)“l  -  l]u/rQ 
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2  Wo 


(3.5) 


Thus,  the  minimum  initial  velocity  is 


°il  ^  (1  ♦  \  Mt)/Mmt  (3.6 

3.2  Forced  Motion  and  Generalized  Forces 
3.2.1  Forced  Motion  Departure  From  m^ 

Consider  that  the  driving  force  p(t)  on  m^  is  directed  toward 
m2,  i.e.,  the  force  vector  is  aligned  along  the  negative  direction  of  vec¬ 
tor,  p.  Therefore,  the  virtual  work  is 


5W  =  j)(  t).6  £_=  -  p(t)6p 

Here,  the  sumbol  "6"  denotes  variation  and  it  should  not  be  confused  with 
the  variable  "6". 

Thus,  the  generalized  forces  in  Equation  (2.8)  are 


Q*  = 

P 


p(t)/m3*u)0  =  -  p*(t) 


0*  =0*  =0*  =0 

vp2  ve  ys 

The  initital  magnitude  of  p(t)  must  be  sufficnent  to  overcome  the  dif 
ference  of  centripetal  and  gravity  forces.  This  gives 

pmtn(0)  =  "3(rro  '  u/rl>  "  V“o(1  ‘  V 
Substitution  of  Equation  (3.1a)  leads  to 

pW01  =  (M2  V'\t  ,3-8 

3.2.2  Forced  Motion  Departure  From  m2 

Again,  the  driving  force  is  considerd  to  be  in  the  direction 
toward  the  end  point,  m^.  The  virtual  work  is 

6W  *  £(t).6Sj  »  p{t)6  |  Pj  -  Pj  |  (3 

Making  use  of  the  constraint  condition  given  by  Equation  (2.1),  yields 

6W  =  p(t)6(i  -  P)  =  -  p(t)6P 

Results  indentical  to  Equation  (3.7)  are  obtained.  The  minimum  magnitude 
of  p(t)  to  get  the  motion  of  m^  started  is 

pm1r>,0)  *  m3(“/r2  ’  r2“o> 

*  "3  t(1  -pM/ro1'2  •  11  -p20/ro)W,'o 


s 


3.2.3  The  Minimum  Energy  Expense  For  Transfer  Operation 

The  first  phase  of  transfer  motion  is  to  drive  across  tht 
system  orbit  and  the  second  phase  is  a  braking  action  which  reduces  velo¬ 
city  of  m^  to  zero  when  It  reaches  the  end.  If  the  velocity  of  at  th< 
instant  it  passes  the  orbit  is  zero,  the  required  energy  has  minimum 
magnitude  which  is  the  sum  of  the  kinetic  energy  given  by  Equation  (3.2) 


and  (3.5).  That  is 


Em1n-2V"'o>2|:(M2t2Mt)2*(lt5V2]/Mmt 

3.2.4  Generalized  Force  Due  To  Air  Drag 


(3.11 


Since  the  transfer  operation  is  accomplished  in  less  than  one 
orbit,  air  drag  has  no  significant  influence  on  the  motion  of  the  system. 
However,  it  can  easily  be  included.  Consider  that  the  drag  force  acts  on 
the  inner  satellite  alone  since  it  is  much  closer  to  earth  atmosphere  if  a 


very  long  tether  is  used.  Let  the  drag  force 


where 


fD  *  cdAv2 


A  =  effective  area  of  m„ 


d  =  air  density  at  altitude  r 2 
c  =  drag  coefficient 


v2  =  Ta  +  T&  +  Tq  (see  Table  1) 


(3.12) 


Then  the  virtual  work  is 


and  results  in 


Q*e  2  -  Pq[Rq  +2(1-  9')]  (3.13) 

★ 

where  =  cdAQ/m.  Note  that  the  air  is  treated  non-rotating. 

3.2.5  Summary  of  Initital  Conditions  and  Generalized  Forces 

To  summarize  the  initital  conditions  of  the  variables  and  the 
generalized  forces  of  all  cases  treated,  Table  3  is  presented. 


TABLE  3  Summary  of  Initial  Conditions  and  Generalized  Forces 


FREE 

Depart  from  m^ 

TRANSFER 

Depart  from  m2 

POWERED 
Depart  from 

TRANSFER 

Depart  from  m2 

6i’  V 

0 

0 

0 

0 

6o 

<1+V 

(1  +  1/2MJ/M 

t  mt 

(1+M3  +  l/2Mt)/Mmt  (1  +  l^Mt)/Mr 

°0 

1  -  Aj/t 

A2/t 

1  -  yi 

A2/t 

‘•iu 

3(1  «• 

none 

none 

x*(0) 

0 

0 

0 

0 

P*,n(0) 

min 

none 

none 

none 

none 

Q*(  t) 

0 

0 

0 

0 

Q*(  t) 

0 

0 

-p*(t) 

-p*(t) 

QJIt) 

-PdCRo+26(1'0,)] 

"PdCr0+26(1"9 ' )] 

-Pd[R0+26 ( 1-0 ' ) 

-Pd[Ro+26(1-9')] 

o;(t) 

0 

0 

0 

0 

3.3  Reversed  Position  of  the  Satellite  System 

Formulation  has  been  based  on  a  configuration  that  the  main 
satellite  is  moving  outside  the  system  orbit.  If  the  two  tether  con¬ 
nected  satellites  Interchange  their  positions,  referred  to  as  the  reversed 
system,  equations  of  motion  can  be  applied  without  changes.  One  simply 


refers  be  the  mass  of  the  outside  satellite.  The  two  systems  have  dif¬ 
ferent  values  for  the  mass  ratios;  the  reversed  system  has  M^>>1  and  M.j<<1 
while  the  regular  system  has  M^<<1  and 
3.4  Initial  Value  of  X* 

To  show  that  A*(0)  =  0,  one  may  use  equations  of  motion  for  zero 

tether  mass  given  by  Fquation  (2.10).  At  an  initital  equilibrium  state: 

9"  =  9  =0,  3"  =9=0,  5*  =o‘  =0,  the  third  and  fourth  equations  are 

o  o 

satisfied  for  any  value  of  A*(0)  =  0.  Now,  setting  the  rest  of  initial 

values,  9*  =  0'  =  0,  a  =1  and  <5  =  (1  +  Md/M  and  \*(0)  =  0  in  the  first 
ooo  o  J 

and  second  equations  and  solving  for  6"  and  a",  yields 

6"  =  3(1  +  M.)/M  and  a"  =  3 
o  3  o 

This  result  satisfies  the  vector  equation, 

^(0)  =  P2(0)  +p  (0) 

=  *(/[-  3(1  +  M,)/M  +  3]i  =  3(M,/M)t«J 

It  has  been  proved  that  the  assumption  \*(0)  =  0  is  correct.  The  inclusion 
of  m.  will  change  the  magnitude  of  5"  and  o"  but  not  the  value  \*(0). 


4.  NUMERICAL  METHODS  OF  SOLUTION 


An  approximate  solution  of  Equation  (2.8),  the  equations  of  motion, 
cannot  be  found  directly  by  any  existing  computer  integration  subroutine 
due  to  the  constraint  condition  on  the  variables.  The  author  has  been 
unsuccessful  in  seeking  a  special  computer  program  for  solution  of  a  system 
of  differential  equations  with  constraints.  A  program  which  combines  an 
integration  subroutine  and  a  subroutine  for  the  determination  of  zeros  of 
analytical  functions  has  been  suggested. 

4.1  Integration  by  Combined  Subroutines 

The  commonly  used  subroutine  "DGEAR"  is  a  differential  equation 
solver  which  finds  approximations  to  the  solution  of  a  system  of  first 
order  ordinary  differential  equations  of  the  form  y^  =  f  (t,  y)  with  ini¬ 
tial  conditions.  The  basic  methods  used  for  the  solution  are  of  implicit 
linear  multistep  type.  The  user  may  use  either  the  implicit  Adams  methods 
(up  to  order  twelve),  or  the  backward  differentiation  formula  methods  (up 
to  order  five),  also  called  Gear's  stiff  methods.  See  Appendix  C. 

4.1.1  Transformation  of  Equations  of  Motion  to  First  Order  System 
To  convert  the  system  of  equations  given  by  Equations  (2.8) 
into  a  first  order  system,  the  variables  are  redefined  as  follows: 


*1 =  6 

y2 

=  yj  =5' 

(4.1a) 

y3  =  0 

y4 

=  y3  =  a’ 

(4.1b) 

y5  =  9 

y6 

-  yj  =  9' 

(4.1c) 

y7  =  13 

y8 

=  y7  =  3’ 

(4. Id) 

Thus,  Equation  (2.8)  can  be  rewritten  in  the  form 


Clly2  +  C12y4  =  Allyl  +  A1 3y 3  +  A16y6  +  A18y8  +  A19X*  +  Qp2  (4,2al 

C21y2  +  C22y4  =  A21y 1  +  A23y3  +  A26y6  +  A28y8  +  A29X*  +  Qp  (4-2bl 

°lly6  +  °12y8  =  B12y2  +  B14y4+  B15y5  +  B17y7+  B19X*  +  Q*9  {4,2cl 

°21y6  +  °22y8  =  B22y2  +  B24y4+  B25y5  +  B27y7+  B29X*  +  (4-2d) 

where  all  the  coefficients  are  defined  in  Appendix  B.  The  constraint 
equation  becomes 

g(y)  =  M2y^  +  (1  +  M3)2y3  -  2M(1  +  M3)y1y3cosy7-  (1  -  y3)  =0  (4.2e) 

Solving  for  y£  ,  y^  from  Equations  (4.2a)  and  (4.2b)  and  yj^  ,  y^  from 
(4.2c)  and  (4. 2d),  results  in 

>’n  ■  S  [  ill  fni^i  *  V*  *  ^  1  n  »  2,  4,  6,  8  (4.3a) 

n 

and  four  more  equations  from  Equation  (4.1), 


y;  »  yn+1  n  =  1,  3,  5,  7  (4.3b) 

where  the  coefficients  are  defined  as  follows: 


=  CUC22  -  Cl2C21,  D6  = 

°8  = 

D11D22  '  D12D21 

(4.3c) 

2i  =  C22A1 i  "  C 1 2  A 2 i 

i 

=  1,  3,  6,  8,  9 

(4.3d) 

4i  =  CllA2i  "  C 2 1 A 1 i 

2i  ’  f4i  "  0 

i 

=  2,  4,  5,  7 

(4.3e) 

6i  =  D  2 2  8 1 i  "  D 1 2  B  2 i 

i 

=  2,  4,  5,  7,  9 

(4.3f ) 

81  =  D 1 1 B  2  i  "  D  2 1 8 1 i 

fcl  =  fOJ  =  0  i  =  1,  3,  6,  8  (4.3g) 


4.1.2  Integration  Procedure  and  Subrouti ne-ZANLYT 


With  the  initial  conditions  yn(0)  for  n  =  1  through  8,  Q*(0) 
and  x*(Q)  given  in  Table  3,  y^x^)  can  be  obtained  by  using  integration 
subrouti ne-DGEAR  provided  that  X*(x^)  is  known.  Let 

xt+i (ii )  =  x T x )  +  axt(Ti)  j  =  1,2..,  k  (a) 


and  rewrite  Equation  (4.3)  in  the  form 


ynj  (tl'  =  fn(,l>  lj(TlH 


Now,  the  above  equation  is  integrated  for  j  =  1,2,:..",  k  by  using  DGEAR  to 
obtain  yn^Ti^****»  ynk^Tl^‘  Substituting  these  in  the  constraint 
equation.  Equation  (4.2e),  yields 

WV'  =  Rj  j  "  1,  2,..,  k  (O 

where  Rj  denotes  the  residual  of  fucntion  g(yn (t ^ ) )  from  zero. 

Subroutine  "ZANLYT"  is  a  program  for  finding  zeroes  of  an 

analytical  function.  It  is  expected  that  ZANLYT  will  determine  X*(x^)  such 

that  as  to  make  the  residual  approximately  equal  to  zero. 

Finally,  X*( x ^ )  is  used  to  integrate  Equation  (4.3)  one  more 

time  to  obtain  yn(x^).  This  procedure  completes  one  integration  step  and 

is  ready  to  trove  forward  to  and  repeat  program  loop. 

There  are  questions  concerning  the  magnitude  of  k,  the  first 
★  ★ 

guess  of  X^x^),  and  magnitude  of  AXj(x^)  to  be  resolved.  Note  that  x*  does 
not  appear  in  Equation  (4.5)  explicitly,  and  nor  there  exists  an  equation 


governing  x*.  More  detailed  study  of  ZANLYT  is  required  to  answer  these 
questions. 

4.2  Elimination  of  the  Lagrange  Multiplier 

Two  possible  approaches  for  the  elimination  of  the  Lagrange 
multiplier  have  been  attemped.  One  is  to  solve  the  constraint  equation  for 


one  of  the  variables  in  terms  of  the  other  three  and  hence  reduce  the 
system  to  three  independent  variables.  Unfortunately,  the  resulting 
expressions  are  so  complex  as  to  make  an  analytical  formulation  of  the 
equations  of  motion  unpractical.  The  other  approach  is  to  eliminate  A* 
directly  from  any  two  equations  of  the  system  of  four  equations;  and  hence 
three  independent  equations  without  A*  can  be  formed.  A  fourth  second 
order  differential  equation  is  obtained  by  differentiation  of  the 
constraint  equation  twice.  Expressing  the  equation  so  obtained  in  a  simi¬ 
lar  form  as  Equation  (4.2),  one  has 


cny2  + 

Cl2y4  =  0 

(4.4a) 

*  C22yi 

'  1  *2iyi  *  ^ab 

(4.4b) 

+  612y8 

"  E  6l(yi  *  ^a 

(4.4c) 

*  °22y8 

'  E  52iyi  *  »cd 

(4.4d) 

where 

5ab  "  A29t)*2  ‘  A19Q; 


°ca  ’  A19%*  -  B19^  •  Myl<V  W! 


P  UU  V  u  ^2 

an'd  the  coefficients  are  defined  in  Appendix  B.  Note  that  the  above 
equations  are  formed  with  the  use  of: 


(4.4a)  =  —2  (4.2e) 
dx 

(4.4b)  =  A2g(4.2a)  -  Aig(4.2b) 
(4.4c)  =  A19(4.2c)  -  Big(4.2a) 
(4.4d)  =  My1(4.2c)  +  M3y3(4.2d) 


Since  second  and  higher  order  terms  were  neglected  from  the  original 


equations  of  motion,  they  are  not  included  in  Equation  (4.4a). 


Now,  Equation  (4.4)  can  readily  be  reduced  to  the  first  order  system, 

n  =  1,  3,  5,  7  (4.5a' 


K  -  Vi 


yn  =  fl J£  fniy1  f  V 


n  =  2,  4,  6,  8 


(4.5b) 


where 


52  =  °4  ~  C11C22  "  C12C21 *  °6  =  °8  =  511°22  "  °12°21 


^2i  =  “  C12S2i‘  f4i  =  CllA2i  1  =  l’  3'  6>  8’  9 


f2i  -  *4i  *  0 


i  =  2,  4,  5,  7 


f6i  =  022Sli  ‘  °12§2i *  f8i  =  5llS2i  '  521Bli  1  =  2>  4*  5’  7*  9 


f6i  -  T8i  -  0 


i  =  5,  8 


^2  =  '  C12Qab’ 


^4  =  CllQab 


^6  =  °22Qca  "  D12Qcd’  ^8  =  5llQca  "  °21Qca 


The  above  system  is  referred  to  as  the  derived  system  while  Equation  (4.3) 
is  called  the  original  system. 

4.3  An  Example  for  Illustration 

An  example  which  has  known  solution  is  given  here  for  illustration 
and  also  may  be  used  for  programming  verification.  A  pendulum  of  mass  m 


and  length  L  with  small  initial  angle  aQ  and  zero  initial  velocity  has  the 


linear  solution  8  =  9qcost  =  ,/gTLt.  If  rectangular  coordinates  x  =  LsinQ 


and  y  =  Lcosa  are  used  for  the  Lagrange  formulation,  one  obtains  the 
equations  of  motion 


with  a  constraint  equation. 
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x  +  y  =  L  or  x*  +  =  0  (4.6c) 

Note  that  Equations  (4.6a)  and  (4.6b)  can  also  be  obtained  from  conditions 
of  equilibrium  as  shown  in  Figure  3(a).  Setting  £  =  x/L,  n  =  y/L,  t  = 
yg7Lt  and  A*  =  A/(mg/L),  the  above  equation  become 

6"  =  A*5  (4.7a) 


>"  =  x*n  +  1 


(4.7b) 


2  2 

5  +  n  =1  or  +  nn'  =  0  (4.7c) 

where  prime  denotes  differentiation  with  respect  to  t. 

To  form  the  derived  system,  one  equation  is  obtained  by  elimi¬ 
nating  A*  from  Equations  (4.7a)  and  (4.7b)  and  a  second  equation  by  dif¬ 
ferentiating  (4.7c).  Thus, 


nC"  -  5n"  =  -  6 


(4.8a) 


«"  - 


nn 


-  (?'  +  n 


■2i 


(4.8b) 


If  a  linearized  pendulum  is  treated,  the  term  on  the  righthand  side  of 
Equation  (4.8b)  can  be  neglected.  Solving  for  5"  and  n"  from  the  above 
equations  results  in 


5"  =-5n  -  US'2  +  n’2) 

2  .2 

n"  =  5  -  nU'  ♦  n  ) 

It  can  be  shown  by  direct  substitution  that  the  following 

S  =  sine  =  9qcost 

.  2 
n  =  cose  =  1  -  1/2(9qcost) 

2  2  2 
A*  =  -  cose  +  (eQsinT)  =  -  [1  -  l£(9  cost)  +  (eosinT)  ] 


(4.9a) 

(4.9b) 


is  an  approximation  to  the  original  system.  Equation  (4.7)  and  the  derived 
system,  Equation  (4.9).  Note  that  A*  is  obtained  from  condition  of 


mg 


(a)  Forces  in  xy-coordi nates  (b)  Forces  in  r  8-coordinates 
Figure  3  Equilibrium  of  Forces 

This  example  serves  the  following  objectives: 

(1)  If  a  computer  program  is  desinged  to  solve  the  original  system, 
the  example  can  be  used  as  a  test  case  to  verify  the  program. 

(2)  This  example  illustrates  that  the  proposed  procedure  for  elimina 
tlon  of  the  Lagrange  multiplier  is  a  valid  approach  for  a  simple 
pendulum. 

4.4  Numerical  Verification 

To  verify  that  the  program  can  find  approximate  solutions  of  a  system 
differential  equations  with  a  constraint  or  that  the  derived  system  can 
truly  represent  the  physical  system,  one  should  carry  out  numerical  com¬ 
putation  of  the  example  given  and  compare  with  the  known  solution.  The 
first  order  systems  of  differential  equations  are  respecti vely: 


(1)  the  original  system 


yj  =  y2  y2  =  X*yl 


y3  =  y4 

y4  =  A*y3  +  1 

g(y) 

2  2 

"  yl  *  *3  ‘  1  '  0 

(4.10) 

deri ved 

system 

yj = 

y2 

y2  -  -  yty3  - 

y3  = 

y4 

yi  ”  yi  -  (y2  *  y4)y3 

(4.11) 

(3)  the  physical  system 

8“  +  sine  =  0  (4.12) 

with  initial  conditions:  8(0)  =  8q  and  8 ' ( 0)  =  0, 

for  small  oscillations.  Equation  (4.12)  is  replaced  by  an  approximately 
nonlinear  differential  equation. 

8*  +  8'  -  C83  *  0  (c  *  1/6)  (4.13) 

which  has  known  approximate  solution, 

1  2  2 

8  =  8qC0S0)T  +  j£C8g( COSu>T  -  cos3u>t  )  +  0(c  )  (4.14) 

where 

»  *  1  -  8eo  9o  <(  1 

Equation  (4.14)  is  sufficient  to  represent  the  exact  solution  of  a  simple 
pendulum  for  small  oscillations. 

Now,  from  solutions  of  Equations  (4.10)  and  (4.11)  compute 

9  (  t  )  =  tan”1^  (t  )/y3(t )] 

and  compare  them  with  Equation  (4.14). 


4.5  Parametric  Study 

The  satellite  system  has  three  mass  parameters,  M^,  and  Mt, 
and  one  length  ratio  Rq  and  their  numerical  ranges  may  be  given  as  follows: 

Mass  of  sub-satellite  /mass  of  main-satellite:  1/20  to  1/5 

Transfer  mass/mass  of  sub-satellite:  1/20  to  1/5 

Mass  of  tether  =  transfer  mass 

Radius  of  satellite  orbit/length  of  tether:  25  to  100 
One  additional  parameter  is  the  magnitude  of  initital  velocity  of  the 
transfer  mass  for  free  transfer  or  the  magnitude  of  the  driving  force  for 
powered  transfer. 

Four  cases  are  to  be  treated:  namely  free  or  powered  transfer  for 
both  Inward  and  outward  transfers.  This  number  is  doubled  if  both  regular 
and  reversed  satellite  systems  are  treated.  If  three  values  (high,  medium, 
and  low)  are  taken  for  each  parameter,  this  gives  8(3*)  =  648  computer  runs. 

4.6  Computer  Print-out  and  Computer  Time  Control 

The  amplitude  of  e(i)  and  e(t)  are  essential  to  the  study  of  sta¬ 
bility  of  the  satellite  system  during  mass  transfer  operations:  u(t)  and 
o'(t)  give  respectively  the  poisiton  and  velocity  of  the  transport  mass. 
These  four  variables  are  required  to  be  ploted  versus  non-dimensional  time 

T  . 

The  integration  variable  for  the  first  order  differential 
equations  is  the  non-dimensional  time  t,  t  =  2 rr  constitutes  one  orbital 
revolution.  One  may  stop  the  computation  when  one  of  the  following  is 
reached: 

(1)  t ^  =  expected  value  for  transfer  operation, 

(2)  |  o'(t)  |  =  e,  where  e  is  some  small  positive  value,  and 

(3)  o(t)  =  Ag/t  when  departure  from  m^  or 

a ( t )  =  1  -  Ai/i,  when  departure  from  m 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

A  dynamical  formulation  of  the  equations  of  motion  of  a  TSS  with  mass 
transport  along  the  tether  has  been  presented.  Four  second  order  differen¬ 
tial  equations  are  obtained  directly  by  the  Lagrange  method  with  one 
constraint  condition  on  the  variables.  A  system  of  four  independent  second 
order  differential  equations  are  derived  from  the  first  by  eliminating  the 
Lagrange  multiplier  among  the  four  equations.  The  fourth  equation  is 
obtained  by  differentiating  the  constraint  equation  twice.  It  is  desired 
to  obtain  computer  solutions  of  both  systems  to  verify  the  validity  of  the 
derived  system. 

There  are  four  cases  for  mass  tranfer  which  are  departures  from  the 
inner  and  outer  sub-satellites  and  free  and  forced  motions  for  each  case. 
Initial  conditions  for  all  cases  are  presented  in  Table  3.  To  cover  whole 
ranges  of  various  combinations  of  parameters,  a  minimum  of  648  computer 
runs  are  required,  and  many  additional  runs  may  be  needed  when  critical 
state  of  equilibrium  arises.  Thus,  systematic  studies  and  records  keeping 
for  such  a  large  volume  data  becomes  a  major  problem. 

The  system  of  equations  of  motion  can  be  reduced  to  very  simple  forms 
if  the  tether  mass  is  disregarded.  It  is  important  to  find  out  how  signi¬ 
ficant  is  the  contribution  of  tether  mass.  No  conclusions  can  be  drawn  on 
how  each  parameter  affects  the  stability  of  the  TSS  until  numerical 
investigations  are  completed. 


Recommendations: 


The  items  (a)  through  (g)  given  in  the  Introduction  may  be  included  in 
a  dynamical  model  one  at  a  time  so  that  contribution  by  each  can  be  eva¬ 
luated.  If  two  or  more  items  are  treated  simultaneously,  analytical  for¬ 
mulation  of  equations  of  motion  becomes  unpractical. 
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A  paperL  J  which  will  appear  in  the  Journal  of  Applied  Mechanics 
suggests  that  a  dynamical  system  with  constraints  may  be  formulated  without 
using  a  Lagrange  multiplier.  This  new  approach  may  provide  another  means  of 


veri fication. 


APPENDIX  A  DEFINITION  OF  NOTATION  USED  IN  TABLES  l  AND  2 


Notations  related  to  kinetic  energy  of  the  tether  in  Table  1  are 
results  of  Equation  (1.12)  and  related  to  potential  energy  of  the  tether  in 
Table  2  are  obtained  from  Equation  (1.14).  Definition  of  these  notations 
are  presented  as  follows: 

Ej  =  2  11  -  M23(l  -  M^)  +  [2  +  M23(l  -  M23)]o} 

H23  -  2M3M23-  2  *  U  *  M3  *  M23  '  2M3M23)o 
E3  .  2  (2  -  3«23(1  -  Mjj)  +  [4  -  3M23(1  -  M23)]o} 

Gj  =  2[1  -  Mj  (1  -  M3)  +  M3U  -  M3)o] 

tJ,“V"23  -  2M3M23 
G4  =  [2  .  M23(l  - 

g5  *  m3u  -  «3)0 

G6  =  [4  -  3M23(1  - 
G,  =  3  [2  -  2M,  (1  -  M,)  +  3M,{1  - 


APPENDIX  B  DEFINITION  OF  COEFFICIENTS  OF  EQUATIONS  OF  MOTION 
The  coefficients  used  in  the  equations  of  (notion.  Equation  (2.8),  are 


defined  as  follows: 


CU  ■  Ma  *  5MtEl 


c12  *  -  Mc  *  5MtE2 


C21  ■  -  <«c  '  6MtE2)/M3 


AU  *  3"a  +  SHt<El  *  e3» 


C22  *  <M3"b  +  KG1  )/M3 


A13  '  3,Mc  *  5MtE2> 


A16  *  '  -  Mc°  +  SHtIEl6  -  E20)1  *18  '  *  2tMc  +  6MtE2>° 

^19  3  ^  "  (1  +  M^oCOSfJ  +  M^.(aa6  -  a  a) 

A21  =  -  13WC  -  gMt(9MM23  -  9(2M  -  1)M3  +  3MRQ  +  3E2  +  G3  +  G4  +  Gg  ]}  /M 

A23  '  3"b  +  KC18M3  *  IG1  *  GS  -  G7)/M3] 

A26  3  -  2(Mb  -  (MMc/M3)6  +  gMt[6M3a2  -  3(2  -  M)Ma6  +  { 3f>IM23 /M3 ) 6 2 ]} 

A28  3  1 2Mb  +  gMt[12M3a  -  3(2M  -  1)6  -  3RQ  +  (2GX  +  2Gg  +  G3)/M3]}a 

A29  =  (M/M3)  {[(1  +  M3)2  -  l]a  +  1  -  M(1  +  M3)6C0SB 

+  ^[M2  62  +  (1  +  MJ)M3a2  -  M(  1  +  2M3)a6  +  2(aba  -  ac6 

°U  ’  V  -  Mc°  *  h(El5  *  E2«>  D12  =  IMc  -  5MtE23a 

D21  *  V  -  (“W4  *  5Mt(Gl°  *  E26)/M3 


°22  S-'V  5WM3)o 


B12  ’  2<"a  *  lMtEl> 


B14  =  ‘  (2Mc  +  5Mt[6M236  ’  3(2  ‘  1/M)M3  ’  3(M  ‘  1)Ro  +  2E2  +  G3a  +  2M 
B15  =  "  3B11  B17  =  G  B19  3  (1  +  M3  +  M^a^asine 


=  -  2  (MM  -  ;M  E  )/M_ 


B„  =  -  3D 


APPENDIX  C  COMPUTER  SUBROUTINE 


IMSL  ROUTINE  NAME:  DGEAR 

PURPOSE:  Differential  Equation  Solver  -  Variable  Order  Adams 
Predictor  Corrector  Method  or  Gears  Method 
USAGE:  Call  DGEAR 
Algori thm 

DGEAR  finds  approximations  to  the  solution  of  a  system  of  first  order 
ordinary  differential  equations  of  the  form  y'  =  f(x,y)  with  initial 
conditions.  The  basic  methods  used  for  the  solution  are  of  implicit 
linear  multistep  type.  There  are  two  classes  of  such  methods  availble 
to  the  user.  The  first  is  the  implicit  Adams  methods  (up  to  order 
twelve),  and  the  second  is  the  backward  differentiation  formular  (BDF) 
methods  (up  to  order  five),  also  called  Gear's  stiff  methods.  In  either 
case  the  implicitness  of  the  basic  formula  required  that  an  algebraic 
system  of  equations  be  solved  at  each  step.  A  variety  of  corrector 
iteration  methods  is  available  for  this. 

DGEAR  and  the  associated  nuclei  are  adaptations  of  a  package  designed 
by  A.C.  Hindmarsh  based  on  C.W.  Gear's  subroutine  DIFSUB. 

ISML  ROUTINE  NAME:  ZANLYT 

PURPOSE:  Zeros  of  An  Analytic  Complex  Function  Using  the  Muller 
Method  With  Deflation 
USAGE:  Call  ZANLYT 
Algorithm 

Muller's  method  with  deflation  is  used. 
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